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Abstract 

We study hydrodynamic oscillations of a cylindrical Bose gas above the Bose- 
Einstein transition temperature using the hydrodynamic equations derived by 
Griffin, Wu and Stringari. This extends recent studies of a cylindrical Bose- 
condensed gas at T = 0. Explicit normal mode solutions are obtained for non- 
propagating solutions. In the classical limit, the sound velocity is shown to be 
the same as a uniform classical gas. We use a variational formulation of the 
hydrodynamic equations to discuss the propagating modes in the degenerate 
Bose-gas limit and show there is little difference from the classical results. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 
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I. INTRODUCTION 



Recently Andrews et al. [|IJ reported a measurement of a sound pulse propagation of Bose- 
condensed cloud in a highly asymmetric cigar-shaped trap. They have measured the sound 
velocity in a Bose-condensate as a function of the density. Zaremba |2| gave a detailed anal- 
ysis of the collective excitations in a cylindrical Bose gas starting from the T = quantum 
hydrodynamic equation of Stringari |3|], which is based on the Thomas- Fermi approxima- 
tion. The sound velocity of a condensate pulse obtained in Ref. is in good agreement 
with the experimental observations More recently, several other theoretical studies have 
discussed the excitations and pulse propagation in a cigar-shaped trap at T = 

In this paper, we consider the analogous modes in a hydrodynamic regime (where col- 
lisions ensure local thermal equilibrium) for a cylindrical Bose gas above the Bose-Einstein 
transition temperature T B ec- Our starting point is the hydrodynamic equations derived by 
Griffin, Wu and Stringari 0. 

II. HYDRODYNAMIC NORMAL MODE EQUATIONS FOR A CYLINDRICAL 

BOSE GAS 

The linearized hydrodynamic equation for the velocity fluctuations v(r, t) derived by 
Griffin, Wu and Stringari is 

m W = 3^7) V(V ' V) ~ V[V ' Wo(r)] " 3 (V ' V ) W °W " ^ Wf/(r ' t} ' (1) 
where Uq(t) is the static cylindrical trap potential and SU(r,t) is a small time-dependent 
external perturbation. The equilibrium local density no(r) and the equilibrium local kinetic 

pressure -Po( r ) m © are given by 

1 k T 

"o(r) = ^3/2(20), P (t) = ^#5/2(>o), (2) 

where £o( r ) = e fcBT ' w is the local equilibrium fugacity, A = (27r/i 2 /ink-sT) 1 / 2 is the 
thermal de Broglie wave-length and g n (z) = Y^iZi zl /l n are the well-known Bose-Einstein 
functions. 



Throughout this paper, we shall limit our discussion to a purely cylindrical harmonic 
trap potential 
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U {t) = ^mu 2 {x 2 + y 2 ). 



(3) 



Inserting this into (|l|), we obtain coupled equations for the radial vi and axial v z velocity 
fluctuations 



<9 2 v_i_ 5Pq ,_ , o— , s 2 ., 

m — - — V_l(V_l • v±) - muj V±(v ± ■ v_l) - -muJ {V ± ■ v±)t ± 



dt 2 3n 



5Po 
3n, 



°T7 ^ 2 \ 

- V_l - -muj r ± 
o o / 



dv z 
dz ' 



(4a) 



m- 



5P 0t7 2 \ dv ± 



(4b) 



dt 2 3n dz 2 \3n / dz 

where we have set 5U (r, t) = since we are interested in normal mode solutions (driven 
solutions will be discussed in Section VI). We use the convention v(r, t) = v w (r)e~ wt . 
Finally, we shall look for solutions of the kind 



v,,, r 



{xf{r ± ),yf(r±),h{r_ L ))e 



ikz 



(5) 



where r± = yfx 2 + y 2 . That is, we assume that the functions / and h do not depend on the 
radial azimuthal angle. 

Using @ in (4), one finds, after some algebra, a coupled set of equations for the radial 
function / and the axial function h 



' d 2 f + _3_ df_\ + ^ 2 /10 + 5^ df_\ _ . k 
dr\ rj_ dr± I I 3 3 dr± I 



c o( r ± 



1 dh 2 ' 

a o U h 

r± arj_ 3 



, (6a) 



oj 2 h = cl(rj_)k 2 h — ik 



The position-dependent local "sound velocity" c (r±) is defined by 

2 5P {r L ) 5k B T 95/2(z ) 
c {r±) = — r = — B(zq), B(z ) 



3mno(r±) 3m 



03/2 (*(>)' 



(6b) 



(7) 
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One can show that the normal mode solutions of (6a) and (6b) satisfy the following orthog- 
onality 

J dr n (r)v*,(r) • v w (r) = 0, if u n > ^ u n , (8a) 

or, more explicitly, 

/ dr ± r ± n (r ± )[rlf:,(r ± )f n (r ± ) + h* n ,(r ± )h n (r ± )\ = 0, if n' ± n. (8b) 
Jo 

Here the label n specifies the different normal mode solutions. We remark that while, it is 
not obvious, if we set the trap frequency ujq to zero, (6a) and (6b) have solutions involving 
Bessel functions J(k±r±) with the expected phonon dispersion relation uj 2 = Cg(/c 2 + k\). 
It is convenient to introduce a dimensionless radial variable 



,. (2k B T\ l/2 



mujn 



(9) 



We also introduce a dimensionless frequency and wavevector 



uj = — , k = kR. 

u 



(10) 



We observe from (2) that R denotes the "size" of the radial density profile produced by 
the harmonic potential trap. In these units, we note z = e /3/10 e~ s and the classical density 
profile n (r±) ~ e~ s . In terms of these new dimensionless variables, it is useful to introduce 
the new functions: 



f(s) = -iRf(r ± ), h(s) = h(r ± ). 



(11) 



Using (6), the coupled equations for / and h are given by 

^/=£i/]-*(f B <*)£-fs), 



u 2 h = -B(z )k 2 h-k 



!b(*)4 + (5b(*)-.)/ 



where we have introduced the operator L 



d 2 f df 



ds 2 



ds 



(12a) 
(12b) 



(13) 



The rest of this paper is based on the equations in (12a) and (12b), which determine the 
normal mode velocity fluctuations using ( |TTD and (|5|). The associated density fluctuations 
8n(r, t) can be found by using the number conservation law 

f = "V'M. (14) 



III. NON-PROPAGATING MODES 

We first examine non-propagating solutions (k = 0) of (12a) and (12b), in which case 
they reduce to the two independent equations 

u 2 f = L[f], iu 2 h = 0. (15) 

There is a trivial zero frequency solution / = 0, h ^ 0, corresponding to = 0, v z ^ 0. In 
this case, the dependence of h on s cannot be determined uniquely. Using fll4f) , one finds 
that 5n = for this mode. The interesting solutions of (|l^) with Co ^ are given by 

h = 0, u 2 f = L[f]. (16) 

These correspond to oscillations only in the radial direction, which we shall now discuss. 

In the classical gas limit, the operator L in ([13]) simplifies since B(zq) = 1. In this case, 
we can obtain the complete set of normal mode solutions of (Jl6|), namely 



Here L n (s) is the Laguerre polynomial defined by 

L n (s) = e s £;(s n e-'), (18) 
and the orthogonal functions are normalized according to (see (8b)) 

d SS e- s f£\s)f^{s) = 5 nn , (19) 
In ordinary variables, the dispersion relation of these k = solutions corresponds to 
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u 2 n = n^ul n = 1,2,3, (20) 

and the associated density fluctuation is given by 

8n(r, t) oc L n (s = r 2 J R 2 ) exp(-r 2 J R 2 )e~ iuJnt . (21) 

In the lowest (n = 1) mode with uj\ = 10u;q/3, fi(s) is independent of s and — 1 — s. 

This is the two-dimensional radial breathing mode. We note that this particular mode 
corresponds to one of the coupled monopole-quadrupole modes found in Ref. |F| for an 
anisotropic trap in the limit that the axial trap frequency (u z ) is set to zero. In this same 
limit (u z = 0), the other mode has zero frequency and a velocity fluctuation is given by 
v w = a(x, y, — 5z). This kind of solution is not described by the form ([5]) which we are 
dealing with in this paper. 

For a degenerate Bose gas, in which B(zq) is now weakly dependent on s through z = 
e^°~ s , one cannot solve (|1(J) analytically. We discuss variational solutions in Section V. 

However one can check that f(s) =constant is a solution, with frequency uj 2 = y^o- This 

F(o) 



is the analogue of the n = 1 normal mode in (|T7|) for the classical gas (we recall that fl is 
independent of s). 

For comparison with (^(]), the analogous non-propagating normal modes in a cigar trap 
at T = have a spectrum given by ]2|U 

^ 2 = 2/(Z + l)^, Z = 0, 1,2, ... (22) 



IV. PROPAGATING MODES 

In this section, we discuss the more interesting propagating solution (k ^ 0) of equations 
(12a) and (12b). In the classical limit (B(zq) = 1), one immediately finds a phonon mode 
solution 
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h = exp(2s/5), / = 0, u 2 = -k 2 . (23) 



This is a longitudinal sound wave with the dispersion relation u = cok, where the sound 
velocity Cq = 5k^T/3m is the same as for a uniform classical gas. There is no radial 
oscillation (i.e. / = 0) associated with this phonon mode in the classical limit. Using (|14|), 
the associated density fluctuation is found to be 



5n(r, t) oc exp(— 3r^/5i? 2 ) exjp(ikz — ic kt). 



(24) 



In the classical limit, one can show that the dispersion relation u = c^k is valid for any 
cylindrical trap potential Uq(t±), i.e. it is not limited to parabolic potentials. In this more 
general case, the solution of the hydrodynamic equation corresponding to the phonon-like 
mode in the classical limit is given by h — exp(2£/o / bk^T) and / = 0, with the associated 
density fluctuation 8n oc exp(— SUo/SksT). 

In order to find how the non-propagating normal mode solutions in (17) are modified 
when k ^ 0, we expand / as follows: 



f(s) = J2aJi°\s). 



(25) 



This follows the approach of Zaremba [0 for a Bose-condensed gas at T = in a cigar-shaped 
trap. Substituting (|25| ) into (12a) and (12b), we obtain the coupled linear equations for the 
coefficients a„: 



u 2 - 



(0) 



2 5 



;k \a n + 



2 k 2 



6 n 



Quj 2 - Ik 2 

where the matrix elements M nn i are defined by 

Mnn ,= f°da s 2 e- s fi°\s)fS\s). 
Jo 

Using the identity for Laguerre polynomials 

/ ds e~ s L n (s)L n ,(s) = 5 nn >(n\) 2 , 
Jo 

we find 



J2 M nn>a n i = 0, 



(26) 



(27) 



(28) 



M nn t = 2nb nn i - \furib n i n ± x - 



(29) 



To lowest order in k 2 , one finds the eigenvalue uo 2 of ( p6|) is given by (M nn = 2n) 

_ 2 10 , /5 M nn \ J2 10 , 13- 2 



k 2 = ^n + ^-k 2 . (30) 



.6 5n J 3 30 

In ordinary units, the excitation spectrum is given by 

= A + ^—k\ n = 1, 2, 3, • ■ • . (31) 

We note that the correction term in (|3lD is of order (kR) 2 relative to the first term, which 
is assumed to be large. Thus this spectrum for propagating modes for a classical gas in a 
cylindrical harmonic trap is only valid for kR <C 1, where R is the radial size of the trapped 
gas density profile. 



V. VARIATIONAL SOLUTIONS 



For a Bose gas above T B ec, one cannot easily solve the coupled equations in (12a) and 
(12b) for k ^ 0. An alternative approach is to recast these equations into a variational form, 



following recent work Jl0| in solving the two-fluid hydrodynamic equations for a trapped 



Bose-condensed gas ||11|| . One finds that the functional 
E[f,h] 



Re j 


rod 

' ds g3/ 2 (zo) 


srL[f] + ^B(z )\h\ 2 k 2 -2sr 


\3 ( " da 3 1 


k 


J 


' o ds g m (z a )(s\f\ 2 + \h\ 2 ) 



• (32) 

has the property that conditions 5E/5f = 0,5E/5h = yields Eqs.(12a) and (12b). Thus 
the normal mode eigenvalues uj 2 are given by the stationary value of this functional E[f, h]. 

For k = with u ^ 0, we can use (^) to estimate the normal mode frequencies using 
the classical solutions of (|T7|) , / = fj® and h = 0, as trial functions. The frequency so 
determined is given by 

ds 93/2(z )sf^L[f^} 



u 2 



jf ds g 3/2 (z )s(f^) 2 

ds [03/2 (*)) - 95/2 



10 10 Jo 

— n H 

3 3 



(,o)]/f(^ + 2^) 



ds g 3/2 (z )s(f^) 2 
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(33) 



The second term in fl33|) gives the quantum correction to the classical limit. 

For k 7^ 0, the most interesting propagating mode is the phonon mode with u oc k. For 
trial functions in fl32|), we take (see (|23|)) 



f = kA f , h = A h exp(2s/5), (34) 

where the constants Af and Ah are real and independent of s. To first order in k, we find a 
phonon-like solution a) = ck, with the (dimensionless) sound velocity c given by 



-2 



g y o c/ 5 / 2 (^o)e 4s/5 — 15 {y ds \9w( z o) - g 5/2 (z )]se 2s/5 J / jf #3/2(^0)5 



g 3/2 (zo)e 4s/5 

(35) 



(36) 



The amplitudes in (33) which are associated with this phonon mode have the ratio 

A f J Q ds l93/2(z ) - g 5/2 (z )}se 2s/5 
h 5 J ds g 3/2 (z )s 

One can see that in a degenerate Bose gas, where (73/2(2:0) 7^ ^5/2 (^0); the radial oscillations 
(Af) are coupled to the axial oscillations (Ah). 

The normal mode frequencies given by the variational expressions in ( |3"3"| ) and ( |33| ) can 
be numerically calculated. All the integrals can be evaluated analytically using the useful 
identity, 

/>nn 

ds g n (z )s m = m\g n+m+1 (z ), (37) 





where z = e^ and g n is the Bose-Einstein function, as defined below (2). It is useful to 
plot the temperature dependence relative to the T B ec f° r an ideal gas in a cigar-shaped trap 
described by (3). For a trap of length L with N atoms, we have (using (2)) 

2ttL r°° 



f ZttL r°° 

N = J dr n (r) = J dr ± r ± g 3/2 (z ) 

nR 2 L 



A 3 Jo 



ds g 3 / 2 (z e s ) 
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(rrujpy/ 2 fk B T\ 



5/2 



05/2 (*oj 



(38) 



When T = Tbec, we have /^o = and hence £0 = 1- Thus (|38f) gives the Bose-Einstein 
transition temperature 



B^BEC 





/2tt^ 2n 


f 1 1 






I C(5/2)_ 



2/5 



(39) 



for a cigar trap in the usual semiclassical approximation. 

In Fig. 1, we show how the sound velocity given by (R3) varies with temperature down to 



TbeC; relative to the classical value c = y5kBT/3m. In Fig. 2, we show the temperature- 
dependent results for the frequencies of the non-propagating modes based on (|3"3"D. We recall 
that since fi is constant, the correction term in (^) vanishes for the n = 1 mode. The 
variational calculations shown in Figs. 1 and 2 indicate that there is little change in the 
normal mode frequencies given by the classical limit for temperature down to T = Tbec- 
As noted at the end of Ref. 0, these results are to be expected since the only place where 
the Bose nature of the gas enters is in the first term of Eq. ([I]). This involves the ratio 
B(zq) = g5/2(zo)/g3/2(zo), which is remarkably close to the classical value of unity, even at 
the center of the trap. 



VI. PROPAGATION OF SOUND PULSES 

In this section, following the approach of Zaremba 0, we discuss the propagation of 
sound pulses induced by a small external perturbation 5U(r,t). We assume that 5U(r,t) 
has no radial dependence and is switched on at t = 0, i.e. the external perturbation is of 
the form 

5U(r,t) = 8U(z)6(t). (40) 

The equation of motion ([[]) with the external perturbation 5U can be solved by introduc- 
ing a Fourier representation of the velocity fluctuations (compare with (^)) and the external 
perturbation 
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v(r,t) 
5U{z) 



—e lkz (xf(k, r ± , t), yf(k, r±, t), h(k, r±, t)), 

-oo ZlT 

^e lkz 6U(k). 
-00 2n 



(41) 



Taking the Fourier transform of ( f4"T| ) and using the radial variable s defined in (|9]), we obtain 
a coupled set of equations (compare with (12a) and (12b)) 
d 2 f 



dt 2 

d ^ + ul{ h -B{ ZQ )W¥-k 



0. 



^ 0)s l + (^ )- s) / 



(42a) 

i—5U(k)S(t). (42b) 
m 



Here we have used notation analogous to that in (0), but now / and h also depend on t 
and are for particular k component. 

In order to solve these coupled equations, we expand / and h as follows, 



f(k, s,t) =J2 b m(k, t)f m (k, s), h(k, s,t) = J2 b m(k, t)h m (k, s) 



(43) 



where the basis functions f m and h m are the normal mode solutions of (12). These have 
eigenvalues O^k) and form an ortho normal set (see (8b)), 



ds g3/2(z )[sf*J n + h* m h n ] = S mn . 
Substituting the expression (f43"D into (42) and using (0), we obtain 

8 2 br 



(44) 



dt 2 



+ uj 2 n (k)b n = F n (k)6(t) 1 



(45) 



where ui n {k) = u n (k)uj and 

F n (k) 



ik $ u ( k ) 1 ds g 3 / 2 (z )h* n (k,s). 



m Jo 



(46) 



The solution of (|4"5| ) with the boundary condition b n (t < 0) = is given by 

F (k) 

b n (k,t)=e(t)^y-sm{u n (k)t], 



(47) 



where the quantum number n refers to the radial degree of freedom. 

In order to analyse the time evolution of the associated density fluctuations, it is conve- 
nient to work with the radially-averaged density 
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5n(z,t) = J dr±8n(r,t). 



(48) 



Using fll4]), one finds 

—5n = - J dr ± V • (n v) = - J dr ± n -^ 



-*-^3 y o ^3/2(^0) 



00 dk 

-00 27T 



ke lkz h(k,r ± ,t) 



nR 2 _ r 00 dfc ikz k 2 SU(k) nu . . ... 



A 3 „ J-00 27r mu) n {k) 
Integrating (f49"D over t, one finds 
vri? 2 



dsg 3 / 2 (z )h n (k, s) 



(49) 



A 3 



E 



00 ^e^^P^(t)[l-cosa; n (A;)t] 



00 



2ty muj^k) 



ds g 3/2 (z )h n (k,s) 



(50) 



One can see that a low frequency phonon mode makes a large contribution to (pOf) . To 
illustrate the contribution to ( pOf) which is associated with the phonon density fluctuations 
uj n = ck in the classical limit, we use h n (k,s) = Ahexp(2s/5) (see (|23f) ) where the normal- 
ization condition (44) gives A\ = 1/5zq. The contribution to (|50|) of this classical sound 
wave is given by (we use ^ = ^^-zo appropriate to the classical limit) 



5n(z, t) = {sU(z) - l -[8U{z - ct) + SU(z + ct)] \ , 



(51) 



which has the form of a propagating pulse moving with a speed ±c. 



VII. CONCLUDING REMARKS 

In this paper, we have given a detailed analysis of the hydrodynamic normal modes of 
a Bose gas in a cigar-shaped trap above T BEC . We have discussed the non-propagating and 
propagating modes, both in the classical limit as well as in the degenerate Bose limit just 
above Tbec- Our results complement the analogous studies @,§-§[ of such modes in the 
quantum hydrodynamic limit at T = 0. In contrast with the T = analysis, which works 
with a single equation for the density fluctuations, we have to work with coupled equations 
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for the velocity fluctuations. For simplicity, we have considered the limit of a uniform gas 
along the axial direction. We note that Stringari [[| has discussed the condensate modes at 
T = in the limit of a very weak trap in the axial direction (u z <C ujq). 

As in Ref. [0, we have ignored the Hartree-Fock mean field contribution to the hydro- 
dynamic equations. Such terms are given in Eq. (6) of Ref. \U\. In place of ([TJ), we obtain 
the linearized velocity equation 

m w = 3^U V(V ■ v) " v[v ■ w(r)] " 3 (v ' v)w(r) 

+2gV( V • n v) - ^8U(r, t). (52) 

This now involves the effective trap potential 

U(r) =U (r) + 2gn (r), (53) 

which also appears in the equilibrium fugacity z = e ^°- u) in the expressions for n (r) 
and -Po( r )- The usual s-wave scattering interaction is g = Anah 2 jvn. The analysis given in 
this paper can be generalized |12| to include the effects of this HF mean field but it is much 



more complicated. We simply quote some final results for the classical limit. The n = 1 
non-propagating mode has a frequency given by 

2 10 2 / gn (r = 0) \ 
"i = (1 - 2kBT ), (54) 

where no(r = 0) is the density at the center of the cylindrical trap. The sound velocity 
corresponding to (23) is given by 

c 2 = SfceT + gn (r = 0)^ 
3m 3m 

The two-fluid hydrodynamic equations for a trapped Bose-condensed gas (T < T B ec) 
have been recently discussed by Zaremba and the authors JO]]. These equations have been 



used to study first and second sound modes in a dilute uniform Bose gas [[L3|] at finite 
temperatures. It is found that first sound corresponds mainly to an oscillation of the non- 
condensate, with a velocity given by 
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2 5 k B T g 5/2 (z ) 2gn 

M, = — r H . (ODj 

3 m #3/2(20) m 

In contrast, the second sound mode mainly corresponds to an oscillation of the condensate, 
with a velocity given by 

2 9 n c0 /rrA 

a, = . 57) 

Here n c0 (n ) is the equilibrium condensate (non-condensate) density. As discussed in 
Ref. [[13], to a good approximation, one can use 

no = -^#3/2(20), (58) 



where the equilibrium fugacity is zq = e 



In principle, we could use the equations in Ref. |TT| to extend the analysis of the present 
paper and discuss the propagating first and second sound modes in a cigar-shaped trap. 
Here we limit ourselves to some qualitative remarks. One expects to find an expression 
similar to (|5"T|) for the propagation of a pulse, and there should be distinct first and second 
sound pulses moving with velocities quite close to u\ and u 2 as given above. However, as the 
expression in ( |5~TD shows, the relative amplitude of these two modes is proportional to 1/uj. 
We conclude that if pulse experiments such as in Ref. were done in the hydrodynamic 
region, most of the weight would be in the second sound pulse if u"l <C u\. This mode, 
given by (|57|) , is the natural hydrodynamic analogue of the Bogoliubov mode exhibited in 
the quantum hydrodynamic region at T = [|],[|. At temperatures close to Tbec, the 
first sound pulse has a much faster speed and thus its intensity will be very weak. The 
observation of distinct first and second sound pulses in cigar-shaped traps would be very 
dramatic evidence for superfluid behavior in dilute Bose gases. The experiment would best 
be done at intermediate or lower temperatures, where u± and u 2 are more comparable in 
magnitude. 
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FIGURE CAPTIONS 

Fig.l: The sound velocity c as a function of temperature relative to Tbec- The values are 
normalized to the classical gas result c = \Jhk-QT/m. 

Fig.2: The normal mode frequencies u n for k = as a function of temperature, as given by 
(|3"3"D. The frequencies are normalized to the radial trap frequency ujq. 
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Figure 1 
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